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I .  INTRODUCTION 


Microwave  interferometry  is  extensively  used  for  the  measurement  of 
travel  and  velocity  of  projectiles  during  their  in-bore  propulsion. 

Though  this  measurement  technique  has  been  with  us  for  a  long  time,  we 
have  not  reached  the  state  where  it  can  be  exploited  for  the  detailed 
analysis  of  projectile  motion  conveniently  and  in  real  time.  In  contrast  to 
projectile  free  flight,  where  microwave  measurement  techniques  are  well 
established  and  Doppler  radar  signals  are  readily  analyzable  with 
stationary  time  series  methods,  the  application  of  microwave  techniques 
to  monitor  projectile  in-bore  motion  is  more  difficult.  This  is  caused 
by  instrumentation  to  operate  in  the  gun  blast  environment  and  by  the 
requirement  for  the  more  complicated  signals  which  present  highly  non¬ 
stationary  time  series  and  contain  a  conglomerate  of  physical  information. 
Various  instrumentation  and  data  analysis  methods  have  been  and  are  being 
developed  to  overcome  the  difficulties.  Basically,  three  types  of  micro- 
wave  instrumentation  are  used  for  the  recording  of  projectile  in-bore 
motion:  the  Michelson  type  microwave  interferometer  with  fixed  frequency, 
the  Doppler  radar,  and  the  microwave  interferometer  with  variable  frequency. 
The  last  one  is  a  recent  development^ . 

Analysis  of  projectile-gun  interface  problems  quite  often  requires 
accurate  knowledge  of  the  projectile  in-bore  and  launch  motion.  The 
reduction  of  projectile-gun  interface  parameters  such  as  pressure,  bore 
traction,  projectile  alignment,  etc.  from  recorded  local  tube  strains, 
is  one  example.  The  analysis  requires  superposition  of  the  elastic  theory 
of  thick  walled  cylinders  on  the  strain  data  using  general  least  squares 
model  fitting  procedures.  One  way  to  reduce  the  dimensions  of  the 
parameter  vector  and  the  constraint  functions  in  the  least  squares  model 
is  the  predetermination  of  the  projectile  motion  by  other  means.  The 
permissible  error  in  projectile  travel  is  less  than  one  millimeter  in 
the  105  mm  M68  tank  gun,  for  instance.  To  determine  whether  microwave 
measurement  techniques  can  be  used  for  such  an  accurate  measurement  of 
projectile  travel,  a  relatively  crude  sensitivity  analysis  was  conducted. 

II.  EXPERIMENTAL  ARRANGEMENT 

The  experimental  setup,  as  customarily  deployed  at  the  Interior 
Ballistics  Division,  represents  a  Michelson  type  microwave  interferometer 
(Figure  2.1).  A  coherent  microwave  beam  polarized  in  the  vertical  plane 
is  radiated  from  the  emitting  antenna  into  the  circular  cylindrical 
waveguide  with  the  free  space  wavelength  A  ,  The  waveguide  propagates 
the  microwave  beam  with  the  guide  wavelength  A  ,  The  beam  strikes 

o 

1H.  Grumann " A  Microwave  Interferometer  With  Variable  Frequency  For 
Continuous  Determination  of  Projectile  Motion ",  Proceedings  of  Workshop 
on  Projectile-uun  Dynamics  under  the  Auspices  of  the  DEA-G-1060  Ballistics 
Research  and  Development ,  24-26  April  1979 ,  USA  ARBADCOM,  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland  21005  (to  be 
publzshed) . 
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Figure  2.1.  Michelson  Type  Microwave  Interferometry  Arrangement  For  The 
Measurement  Of  Projectile  In-Bore  And  Launch  Motion 


the  moving  boundary  and  is  reflected  backward  through  the  waveguide  and 
free  space  to  the  receiving  antenna.  There,  the  incident  wave  is  brought 
to  interference  with  the  outgoing  wave  and  the  resulting  phase  shift  is 
recorded.  The  argument  of  the  sinusoidal  function  varies  between  zero 
and  2tt,  depending  on  the  momentary  superposition  of  the  interfering 
waves.  As  long  as  the  optical  path  length  of  the  microwave  stays 
constant,  the  phase  shift  also  remains  constant.  When  the  path  length 
changes,  a  phase  difference  occurs.  In  particular,  when  the  path  length 
of  the  beam  changes  by  one  wavelength,  the  phase  will  shift  by  2tt.  Because 
of  the  geometrical  arrangement,  a  change  in  path  length  by  one  wavelength 
corresponds  to  a  displacement  of  the  moving  boundary  by  a  half  a  wave¬ 
length,  thus  allowing  us  to  relate  the  change  in  phase  to  the  axial 
motion  of  the  boundary. 

III.  SENSITIVITY  ANALYSIS,  GENERAL 

In  the  experimental  arrangement,  the  waveguide  is  the  gun  tube  and 
the  moving  boundary  is  the  front  of  the  projectile.  The  wavelength  for 
an  infinite  circular  cylindrical  waveguide  is  given  by^ 


2 

D.E.  Gray,  American  Institute  of  Physios  Handbook ,  3rd  ed.  ,  McGraw- 
Hill  Book  Company,  N.Y.  1972. 
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(3.1) 


Ag  =  X0/Il  -  Uo/A*)2]2,  with  Xq  =  c/(nf)  and  A*  =  2irap/R, 

where  Ag  is  the  guide  wavelength,  A0  is  the  free  space  wavelength.  A*  is 
the  cutoff  wavelength  related  to  free  space,  c  =  (y0e0)-*  is  the  velocity 
of  light  in  vacuum,  f  is  the  continuous  wave  frequency,  n  =  (p'e')?  is 
the  refractive  index,  e'eQ  is  the  permittivity,  is  the  permeability, 

a  is  the  radius  of  the  waveguide,  p  is  a  correction  factor  to  account  for 
the  helical  bore  configuration  introduced  by  the  rifling,  and  R  =  jn  ^  or 
in,£  is  the  &-th  zero  of  the  Bessel  function  Jn(z)  or  its  derivative’ Jn(z) , 
respectively,  depending  on  the  propagating  waveguide  mode. 

These  parameters  may  not  stay  constant  during  the  interior  ballistic 
propulsion  cycle.  To  analyze  the  effects  of  parameter  variation  on  the 
guide  wavelength,  we  may  carry  out  an  error  budget  analysis.  Variation 
of  the  parameters  in  the  equation  for  the  guide  wavelength  yields 

AAg/Ag  =  -  [1/(1  -  S2)][Af/f  +  An/n] 

-  [S2/  (1  -  S2)][Aa/a  +  Ap/p  -  AR/R] , 


S  =  (Vy  •  (3.2) 


Since  reflection  of  the  microwave  on  the  projectile  front  may  excite 
higher  propagating  waveguide  modes,  R  must  be  considered  as  a  variable. 

This  equation  shows  that  uncertainties  in  the  frequency  and  the  refractive 
index  are  increased  by  [1/(1  -  S2)J,  whereas  uncertainties  in  the 
parameters  of  the  cutoff  wavelength  are  amplified  by  JS2/(1  -  S2)]  .  The 
functional  behavior  of  these  two  amplification  factors  is  shown  in  Figure 
3.1.  To  illustrate  the  influence  of  the  microwave  frequency  selected  on 
the  amplification  of  errors,  values  for  S  characteristic  to  two  in-bore 
microwave  interferometry  techniques  are  superimposed.  The  points  marked 
by  TE*11  and  TM*j  on  the  g(s)  curves  correspond  to  a  technique  in  which 
the  microwave  frequency  is  chosen  such  that  only  one  propagation  mode  can 
be  excited  in  the  gun  tube.  The  lower  and  upper  frequency  bounds  given 
are  the  cutoff  frequencies  for  the  TEU  and  TE21  modes  and  the  TMol  and 
"^11  modes ,  respectively.  Since  we  are  at  the  present,  mainly  concerned 
with  the  105  mm  tank  gun  system,  the  sensitivity  analysis  will  be 
arbitrarily  restricted  to  this  caliber,  thus  fixing  the  waveguide  radius. 
From  Figure  3.1  we  can  see  that  variations  in  the  frequency  or  refractivity 
and  in  the  parameters  of  the  cutoff  wavelength  are  propagated  into  the 
guide  wavelength  with  amplifications  of  ^  1.5  and  £  0.6,  respectively,  for 
the  single  mode  technique.  The  values  marked  15  GHz,  105  mm  on  the  g(s) 
curves  refer  to  the  in-bore  microwave  interferometry  arrangement  shown 
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Figure  3.1.  Amplification  Of  Errors  In  The  Guide  Wavelength  Due 
To  Variations  Of  Parameters  In  The  Equation  For  The 
Guide  Wavelength 

in  Figure  2.1,  employing  a  15  GHz  microwave  beam  which  is  linearly 
polarized  in  the  vertical  plane.  By  going  to  a  higher  frequency,  the 
amplification  factor  can  be  reduced  substantially.  Moreover,  an  increase 
in  the  carrier  frequency  well  beyond  the  cutoff  frequency  of  the  next 
higher  propagation  mode  results  in  the  excitation  of  a  multitude  of 
propagating  modes  in  the  waveguide.  Depending  on  the  reflecting  surface 
most  of  the  available  energy  will  be  partitioned  into  the  lowest  mode, 
thus  providing  us  with  an  apparent  wavelength  which  fluctuates  between 
two  bounds.  Should  we  be  able  to  design  an  analysis  method  which  permits 
us  to  extract  the  contribution  from  the  lowest  mode  from  the  microwave 
interferometer  data,  the  higher  frequency  approach  would  be  desireable 
for  minimizing  the  error. 

To  obtain  realistic  estimates  for  the  variation  of  the  parameters 
occurring  in  the  equation  for  the  guide  wavelength,  two  different  sets 
of  data  will  be  used.  They  will  be  referred  to  as  shell  push  data  and 
M68  firing  experiment  data. 
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During  1973  a  few  microwave  interferometry  measurements  were 
carried  out  at  the  Interior  Ballistics  Division  of  the  Ballistic  Research 
Laboratory  by  J.  W.  Evans,  where  projectiles  with  different  nose 
configurations  were  pushed  through  a  shortened  105  mm  rifled  M68  cannon 
with  a  shell  pusher  machine.  Three  microwave  frequencies,  10  GHz,  35  GHz 
and  70.65  GHz  were  used.  The  interference  patterns  for  the  carrier 
frequencies  were  recorded  as  functions  of  projectile  displacement  on  a 
chart  recorder.  In  the  analysis,  data  from  test  #3  (11  Sept.  1973)  which 
employed  a  flat  nosed  projectile  and  10  GHz  and  70.65  GHz  microwave 
simultaneously  are  used. 

The  second  set  of  data  was  collected  during  Fall  and  Winter  of  1977 
by  a  well  instrumented  projectile-gun  firing  experiment3  which  addressed 
the  interior  ballistic  causes  of  the  erratic  flight  behavior  of  M392  type 
of  projectiles  when  fired  from  a  105  mm  M68  tank  gun. 

IV.  VARIATION  IN  FREQUENCY 

It  is  common  practice  to  use  the  frequency  as  it  appears  on  the 
frequency  meter  of  the  microwave  generator  as  the  real  microwave  frequency. 
However,  this  may  not  be  true  because  the  frequency  indicator  may  not  be 
exactly  calibrated  to  the  frequency  of  the  microwave  generator  or  the 
frequency  may  slightly  shift  during  the  operation  of  the  instrumentation. 
With  today's  technology,  we  can  safely  assume  that  the  last  digit  in  the 
dial  setting  is  correct  within  ±  0.5,  regardless  of  the  frequency  selected, 
for  our  15  GHz  radar,  that  translates  into  an  uncertainty  in  frequency  of 
I  Af/f|  £  0.00033.  The  resulting  error  jn  the  wavelength  for  the  TBii 
mode  is  about  the  same,  | AA  /A  |  £0.00034.  Hence,  the  maximum  uncertainty 
in  projectile  travel  due  to  this  uncertainty  is  on  the  order  of  0.7  mm 
for  every  meter  in  projectile  travel.  Generally,  the  error  should  be 
smaller . 

The  uncertainty  in  the  carrier  frequency  can  completely  be  eliminated 
by  monitoring  the  microwave  frequency  with  a  counter  throughout  the 
experiment  -  a  procedure  which  will  also  take  out  human  errors  in 
correctly  setting  or  reading  the  frequency  dial. 

V.  VARIATION  IN  REFRACTIVITY 

The  variability  of  the  refractivity,  N  =  (n  -  1)  106,  for  microwaves 
in  air  has  been  thoroughly  investigated.  It  can  be  described  by4 


R.K.  Loder  and  J.O.  Pilcher,  ", Nondestructive  Test  Method  To  Establish  the 
Performance  of  Projectile-Gun  Systems, "  Proceedings  of  the  26-th  Defense 
Conference  On  Nondestructive  Testing,  J5  to  17  November  1977,  Seattle, 
Washington;  Published  by  the  US  Army  Materials  and  Mechanics  Research 
Center,  Watertown,  Massachusetts. 

f 

B.R.  Beau  and  E.J.  Dutton,  Radio  Me teoro logy ,  National  Bureau  of  Standards 
Monograph  92,  1  March  1966. 
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(5.1) 


N  =  (n  -  1)  106  =  KjPp/T  +  K2Pw/T  +  ^/T2, 


where  and  K3  are  constants,  T  is  the  temperature  in  degrees  Kelvin, 

and  Pp  and  P^  are  the  partial  pressures  of  dry  air  and  water  vapor  in  mbar, 
respectively.  Using  the  values  for  the  constants  as  recommended  by 
Smith  and  Weintraub  ,  the  above  equation  can  be  expressed  by 


N  =  (77.6  P  -  5.6  Pw  +  3.75  *  105  Pw/T)/T,  (5.2) 


where  P  is  the  total  air  pressure,  P  =  PD  +  P^.  As  long  as  one  stays 
away  from  frequencies  near  the  water  vapor  resonance  (f  =  22.235  GHz, 

A^  =  13.5  mm)  and  the  oxygen  resonance  (f  =  60.000  GHz,  A0  "  5  mm),  the 
dispersion  of  the  refractive  index  is  within  the  .5%  accuracy  limit  of 
the  above  equation.  This  equation  for  the  variability  of  the  refractivity 
has  been  derived  for  atmospheric  pressures  and  temperatures  up  to  250°C. 

Its  applicability  for  estimating  the  refractive  index  of  the  compressed 
air  column  in  front  of  the  projectile  is  highly  questionable,  even  in 
the  absence  of  blowby  gases.  However,  due  to  the  lack  of  a  suitable 
approximation  or  data,  the  above  equation  will  be  used  for  estimating 
the  order  of  magnitude  in  the  variation  of  the  wavelength  due  to  the 
variation  in  the  refractivity,  utilizing  data  obtained  from  the  M68 
firing  experiment. 

With  the  data  given  in  Table  5.1,  we  can  estimate  the  variability 
of  the  refractive  index,  using  equation  (5.2),  and,  consequently,  the 
change  in  wavelength  for  the  TEji  wave  mode  propagating  in  the  compressed 
air  column  (Table  5.2).  The  uncertainty  in  projectile  travel  due  to  the 
change  of  refractivity  in  the  compressed  air  column  is  on  the  order  of  a 
few  millimeters  and  can  normally  be  neglected.  However,  it  exceeds  the 
less  than  one  millimeter  accuracy  criterion  for  projectile  travel,  that 
is  required  for  the  analysis  of  local  tube  strains.  The  estimate  of  the 
change  in  refractivity  due  to  the  compressed  gas  in  front  of  the  projectile 
neglected  chemical  reactions  between  the  molecular  species,  electron 
excitation,  activation,  dissociation,  etc.  and  did  not  include  blowby 
gas.  The  actual  error,  therefore,  could  be  considerably  larger  than 
estimated  here. 

Figure  5.1  depicts  wavelength  versus  position.  The  uppermost  line 
corresponds  to  a  guide  wavelength  with  the  refractive  index  set  equal  to 
one,  which  is  customarily  done  in  interior  ballistics.  This  vacuum 
approximation  introduces  an  error  in  the  guide  wavelength  and  projectile 
travel  on  the  order  of 


^E.K.  Smith  and  S .  Weintraub n The  Constants  In  The  Equation  For  Atmospheric 
Refractive  Index  In  Radio  Frequencies" 3  Proc .  IRE413  p  1035-1037. 
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TABLE  5.1.  DATA  PERTAINING  TO  THE  ANALYSIS  OF  VARIATION  IN  REFRACTIVITY 


ATMOSPHERIC  DATA: 

Tq  =  300°K,  Pq  =  1.058  *  10^  Pa,  P  =  3.8  *  10^  Pa  (corresponding 
to  90%  relative  humidity) 

MICROWAVE  FREQUENCIES: 

fj  =  15  GHz  and  f2  =  2.725  GHz 

WAVEGUIDE  PARAMETERS: 

2a  =  105  mm,  p  =  1.015,  R  =  1.84  (corresponding  to  the  TE^  mode) 
SHOCK  FRONT  AND  PROJECTILE  MOTION  DATA,  AVERAGED  OVER  A  POPULATION  OF 


FIVE  ROUNDS 

AND  COMPUTED 

FROM  PRESSURE 

DATA  FOR  A  PROOF 

SLUG: 

1) 

Xs,  [m] 

2) 

AXsp,  [m] 

3) 

C 

 -~5y- 

7T  T 

- 6T“ 

TT 

S 

r> 

M 

s 

3.050 

.310 

7.23 

20.8  2.88 

19.4 

2.88 

4.130 

.495 

6.69 

27.8  4.16 

25.2 

3.45 

5.223 

.650 

6.67 

34.8  5.14 

30.7 

3.91 

1) 


2) 


X  .  . 

s 


AX 


sp 


.location  of  shock  front  in  cannon,  corresponding  to  pressure 
gauge  positions  12,  16  and  20,  respectively;  distance  is  measured 
from  breech  end  of  cannon 

length  of  compressed  air  column  in  front  of  projectile 


D/Do...mean  density  of  compressed  air  column, 


4^7  = 
s>7  - 


(X  -  X  )/AX  = 
s  po'  sp 

dimensionless;  X  =  .821  m  is  the  rest  position  of  the  projectile 

nose  in  gun  tube;  D  is  the  mean  density  of  the  compressed  air 
column;  and  Dq  is  the  density  of  the  ambient  air. 

P/Po...mean  pressure  of  compressed  air  column,  dimensionless. 


6) 

7) 


TT  = 
S 


tt/£  =  T/To...mean  temperature  of  compressed  air  column,  dimension¬ 
less;  humid  air  is  assumed  as  an  ideal  diatronic  gas. 

P  /P  ...measured  pressure  immediately  behind  the  shockfront, 
dimensionless . 


Ms...Mach  number  for  shock  velocity 
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TABLE  5.2.  VARIATION  OF  GUIDE  WAVELENGTH  AND  PROJECTILE  TRAVEL  DUE  TO 
THE  CHANGE  OF  THE  REFRACTIVITY  IN  THE  COMPRESSED  AIR  COLUMN 


X 

[m] 

N  *  106 

An/no 

[%] 

AA  ./A  , 

gl  gl 

m 

4V/Jg2 

[’*] 

Ax^ 

[mm] 

Ax2 

[mm] 

3.050 

2300 

.  188 

-.190 

-.296 

1.18 

1.84 

4.130 

2015 

.  159 

-.161 

-.251 

1.60 

2.49 

5.223 

1961 

.154 

-.156 

-.243 

2.03 

3.  15 

Figure  5.1.  Variation  Of  Guide  Wavelength  Due  To  Variation  In  Refractivity 
For  The  105  mm  M392  Projectile-M68  Cannon  System;  Microwave 
Parameters:  15  GHz,  TE^  Propagating  Mode 
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0.085%  ,  and 


AA  ,/X  ,  “  +  0.043%  ,  Ax./x..  =  - 

gl  gl  11 


-  +  0.066%  ,  A*2/x2  ~  -  0.133% 


g2  g2 


(5.3) 


for  the  two  carrier  frequencies  f^  and  f2,  respectively.  The  accumulated 
errors  in  travel  over  the  total  length  are  about  4  mm  and  6.7  mm, 

«  respectively.  We  can  eliminate  this  type  of  error  by  monitoring  atmos¬ 

pheric  pressure,  temperature  and  humidity  during  the  experiment  and 
accounting  for  them  in  the  data  analysis  or  by  measuring  the  refractivity 
.  of  the  ambient  air  with  a  refractometer .  The  next  lower  line  in  Figures 

5.1  corresponds  to  the  guide  wavelength  with  the  refractive  index  adjusted 
to  the  atmospheric  environment.  It  represents  an  upper  bound  for  the 
actual  guide  wavelength.  Estimated  mean  values  for  the  three  selected 
shock  positions  are  shown  underneath.  At  muzzle  exit,  the  guide  wave¬ 
length  converts  back  to  the  free  space  wavelength.  Its  functional  behavior 
can  be  computed  from  the  theory  of  finite  wave  guides,  in  principle. 

Because  of  the  shock  formation  and  the  highly  compressed  air  in  front  of 
the  projectile,  uncertainties  in  the  wavelength  will  remain  for  the 
transitional  ballistic  regime. 

VI.  VARIATION  IN  CUTOFF  WAVELENGTH 

The  variation  of  the  guide  wavelength  due  to  the  variation  of  the 
cutoff  wavelength  is  given  by 

AAg/Ag  =  -  [S2/  (1  -  S2)][.Aa/a  +  Ap/p  -  AR/R]  .  (6.0.1) 


6.1  Variation  of  Bore  Diameter 


The  variation  of  the  bore  radius  is  caused  by  tube  wear.  Taking 
the  cannon  wear  profile  as  obtained  in  the  M68  firing  experiment  (Figure 
6.1),  we  can  approximate  the  variation  in  the  bore  radius  as  function  of 
projectile  nose  displacement  from  the  rest  position  by 


0.0  m  <_  Ax  £  1.0  m:  Aa/a  »  0.0085  (1  -  Ax) 


1.0m  <  Ax 


Aa / a  =  0.0  , 


(6.1.1) 


and  consequently,  estimate  the  variation  in  the  guide  wavelength  to  he 

AY'g  =  "  °-°°85  [S2/(1  "  S2)](1  '  Ax)-  (6,1, 2) 
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WEAR  ten) 


Figure  6.1.  Diametral  Wear  Of  The  105  mm  M68  Tank  Cannon,  No.  10328. 


The  resulting  error  in  projectile  travel  after  one  meter  of  in-bore 
travel  is  0.1  mm  and  4.9  mm  for  the  two  frequencies  f ^  and  f„,  respectively. 

The  105  mm  M68  tank  cannon  example  shows  that  bore  wear  should  be  accounted 
for  in  large  caliber  guns,  if  we  employ  the  single  mode  technique.  The 
multimode  approach,  on  the  other  hand,  allows  us  to  minimize  the  error 
introduced  by  bore  wear  by  increasing  the  carrier  frequency  sufficiently. 

6.2  Effect  of  Bore  Rifling 

Because  most  of  our  gun  tubes  have  a  rifled  bore,  we  must  expect 
that  the  guide  wavelength  equation  derived  for  a  smooth  circular 
cylindrical  waveguide  is  only  a  first  order  approximation.  A  rigorous 
treatment  must  account  for  the  helix  type  of  boundary  conditions  introduced  ' 

by  the  rifling.  Experiments  conducted  to  determine  the  effect  of  para¬ 
meter  variations  on  the  equation  for  the  guide  wavelength  suggest  that 
the  geometric  guide  radius  has  to  be  replaced  by  an  "electromagnetic 
radius",  thus  changing  the  cutoff  frequency  by  a  factor6.  The  correction 


L.  Fortun ,  "Snelheidsmeting  aan  projektielen  in  en  direkt  buiten  een 
kanon  met  behulp  van  Doppler  radar/'  Handleiding  Commission  van  Proefneminq 
deel  XAj  hoofdstuk  5C-1 ,  s  'Graverihage,  July  1963. 
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7 

factor,  p,  must  be  determined  experimentally.  Dr.  Biele  suggested  an 
electromagnetic  diameter  which  is  "approximately  given  by  the  nominal 
caliber  of  the  gun  plus  1.5  times  the  groove  depth"  for  the  TMQ1  and  TE1X 
modes.  This  formula  can  be  translated  into  an  equation  for  p. 


p  =  1  +  k  (6a/2a) 


(6.2.1) 


where  k  is  the  multiplication  factor  and  6a  is  the  groove  depth.  This 
formula  yields  a  value  of  p  =  1.0163  for  the  105  mm  tank  gun.  Experimental 
electromagnetic  diameters  obtained  by  Tengeler8  for  the  TE.  mode  cover 
this  value  (Table  6.1). 


TABLE  6.1.  CORRECTION  FACTOR  BASED  ON  TENGELER'S  EXPERIMENTAL  VALUES  FOR  ELECTROMAGNETIC  DIAMETERS 


CANNON 

DIAMETER  (mm) 

ELECTROMAGNETIC 
DIAMETER  (mm), 
TEn  MODE 

CORRECTION 
FACTOR  p 

MULTIPLICATION 
FACTOR  k 

VALUES  APPLIED  TO 

105  mm  M68  TANK  GUN 

RIFLING 

BORE 

MIN.  p 

MAX .  p 

40L60  NEW 

41 . 10 

40.03 

40. 83± . 15 

1 .0200*. 0037 

1 .50* . 28 

1.0133 

1.0194 

WORN 

(40.84) 

40.89t.15 

1.0215*. 0037 

1.61*. 28 

1.0145 

1.0206 

25  pdr  NEW 

89.6 

87.6 

88.75*. 15 

1 .0131*. 0017 

1 .15*. 15 

1.0109 

1.0142 

WORN 

(88.4) 

00 

1  + 

NJ 

1 .0171*. 0023 

1.501.20 

1.0142 

1 .0185 

To  determine  a  more  accurate  value  for  the  105  mm  M68  tank  gun  ,  an 
attempt  was  made  to  estimate  the  electromagnetic  diameter  from  the  105  mm 
shell  push  data.  Unfortunately,  the  microwave  interferometry  shell  push 
experiment  was  only  set  up  for  the  purpose  of  measuring  the  in-bore  radar 
cross  section  of  preselected  projectile  nose  configurations  and  not  for 
a  detailed  in-bore  microwave  propagation  study.  Since  the  shell  pusher 
machine  had  been  broken  and  was  not  available  for  conducting  the  required 
measurements  with  all  parameters  recorded,  the  data  were  analyzed  to 
obtain  the  desired  information. 

The  data  were  recorded  on  chart  paper;  thus,  they  eluded  application 
of  modern  data  analysis  techniques  such  as  numerical  filtering,  spectral 
analysis,  and  general  least  squares  model  fitting  for  determining  the 
guide  wavelengths  and  amplitudes  of  the  individual  propagating  modes. 
Because  of  the  large  error  involved  with  manually  taking  values  from 
charts,  data  averaging  was  resorted  to. 

-  - 

J.K.  Btele 3  "Measurement  of  In-bore  Motion  of  Projectiles  And  Simultaneous 
Data  Transmission  From  Built-In  Sensors  By  Means  of  Microwave  Inter¬ 
ferometry 3  "  Proceedings  of  the  2nd  International  Symposium  on  Ballistics 
9-11  March  1976 ,  Daytona  Beach ,  Florida. 

Q 

T.I.  Tengeler } "  Het  bepalen  van  de  elektrische  diameter  van  kanonnen3  " 
Leok-rapporten  No.  14.099  en  14.258 3  Oegstgeest,  March/April  1963. 
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Since  the  experiment  was  carried  out  in  an  air  conditioned  room  the 

environment  was  assumed  as  1015  mb  pressure  (sea  level},  25°C  temperature 

and  60%  relative  humidity,  yielding  a  refractive  index  of  n  =  1.000366. 

o 

The  frequencies  of  the  two  carrier  microwaves  as  read  from  the 
frequency  meters  of  the  free  running  oscillator  were  assumed  to  be  correct. 
Taking  only  those  data  for  which  the  projectile  nose  protruded  from  the 
muzzle,  the  following  values  for  the  two  nominal  frequencies  were  obtained: 


f1  =  70.65  GHz:  ^  =  4.2343  mm,  ^  =  c/CfjA  )  =  1.00214 
^  =  c/fn^)  =  70.775  GHz 


f2  =  10.00  GHz:  \2  =  29.9430  mm,  n2  =  c/(f2X2)  =  1.001211 

f2  =  c/(nor2)  =  10.0084  GHz  .  (6.2.2) 


Though  the  frequencies  are  within  the  accuracy  of  the  measuring  device, 
the  observed  discrepancies  between  the  refractivities  and  between  the 
frequencies  suggested  that  the  frequencies  or  the  length  scale  or  both 
are  not  quite  accurate.  Since  the  measurements  were  made  concurrently, 
adjustments  can  be  made  by  minimizing  the  function 

F  =  (fj_  “  ^)2  +  (f2  3?2)2,  £  -  c/(nQkT)  .  (6.2.3) 


The  length  scaling  factor  k  is  then  determined  by 

k  =  (c/nQ)  £ (1/X1  )2  +  (l/^)2]/!^/*!  +  f2/X2]  =  1.00175  ,  (6.2.4) 


yielding  the  following  values  for  the  free  space  wavelengths  and 
frequencies  for  the  two  carrier  microwaves: 


X^  =  4.2417  mm 


X2  =  29.9954  mm 


^  =  70.6514  GHz. 


^2  =  9.9910  GHz  . 


(6.2.5) 
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The  resulting  adjustments  in  frequency  were  Q.002%  and  0.090%,  respectively. 

I  n  order  to  reduce  the  effect  of  higher  propagating  modes  on  the 
wavelength  for  the  most  dominant  mode,  the  dominant  mode  was  determined 
as  well  as  the  length  of  the  interval  required  for  data  averaging.  For 
that  purpose,  the  wavelengths  and  interference  cycles  for  the  first  few 
propagating  modes  were  estimated  from  the  already  known  information 
(Tables  6.2  and  6.3),  setting  p  =  1.015.  The  interference  cycle  gives 
the  number  of  sinusoidal  cycles  necessary  to  produce  a  path  change  of 
one  wavelength  between  the  two  waves. 

The  interference  cycles  marked  with  •  and  ■  refer  to  the  TE;q  and 
TMqi  mode,  respectively,  as  the  reference  mode.  The  experimental  data 
of'  the  shell  push  experiment  cover  approximately  600  and  100  cycles  for 
the  two  frequencies  f^  and  and  are  given  in  Figures  6.2  and  6.3, 
respectively. 

TABLE  6.2.  INTERFERENCE  CYCLES  PRODUCED  BY  PROPAGATING  MODES  FOR  THE  CARRIER  FREQUENCY 

n, 

fj  =  70.6514  GHz 


Xg.  [mm] 

MODES 

TEU 

TM  . 
ol 

TE21 

TE  .  TM. . 
ol  11 

TE31 

4.24287 

TEn 

x 

4.24378 

™ol 

5199 

4.24487 

te2i 

2116  • 

3568 

4.24672 

TEol 

™11 

1101  • 

1397  ■ 

2298 

x 

4.24773 

TE31 

872  • 

1048 

1484 

4192 

X 

4.25072 

™21 

540 

603  ■ 

726 

1062 

1422 

TABLE  6.3.  1 NTE REERENCK  CYCLES  PRODUCED  BY  PROPAGATING  MODES  POR  THE  CARRIER  FREQUENCY 

% 

f2  =9.9910  GHz 


o 


(““)  HX9N3T3AWM 
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Figure  6.2.  Wavelength  Versus  Projectile  In-Bore  Displacement  For  A  Microwave  Frequency  Of  70.6514 


32,244  h-  6UI0E  WAVELENGTH  (TEn MODE) 
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Figure  6.3.  Wavelength  Versus  Projectile  In-Bore  Displacement  For  A  Microwave  Frequency  Of  9.9910 


The  modulation  of  wavelength  and  amplitudes  indicates  that  about  one 
third  of  the  microwave  energy  is  partitioned  into  the  higher  propagating 
modes.  To  remove  this  modulation,  proper  averaging  of  the  data  is 
required.  This  is  not  possible  for  the  first  data  set,  because  the 
interference  cycles  produced  by  the  higher  modes  exceed  by  far  the 
number  of  periods  available  for  averaging.  However,  for  the  second  data 
set  C±2  =  9*991  GHz)  the  situation  is  reversed.  I  f  we  select  99  periods 
as  the  length  for  our  data,  the  interval  used  for  averaging  will  be 
factors  of  2.5,  5.0,  and  6.5  larger  than  the  periods  for  the  TE ^  ~  TE^, 
TE  ^  -  TE-j^,  and  TE^  -  TE^  mode  interferences,  respectively. 

Averaging  the  data  over  this  interval,  we  obtain  for  the  guide  wavelength 


A  =  30.41414  ±  0.00192  mm  . 


(6.2.6) 


The  error  given  corresponds  to  one  standard  deviation.  Equation  3.1  can 
be  rewritten  as 


p  =  [R/(2ira)J[c/Cn*)]/{l  -  [c/ (n^A)] 2 }1/2  ,  (6.2.7) 


expressing  the  correction  factor,  p,  as  function  of  parameters  which  have 
been  determined  already.  Assuming  that  p  is  constant  over  the  entire 
tube  length  and  inserting  the  appropriate  values  for  the  parameters 


R(TEn)  =  1.84118,  2a  =  0.105  m  , 

n  =  1.000366,  £  =  9.991  *  109  Hz, 

we  obtain  the  following  estimate  for  the 


c  =  2.997925  *  IQ8  m/sec, 

and  T  =  0.03041414  m  ,  (6.2.8) 

correction  factor 


p  =  1.01227  . 


(6.2.9) 


This  value  is  in  excellent  agreement  with  p  =  1.01255  ±  0.00165  obtained 
from  Tengeler’s  data  for  a  new  87.6  mm  cannon,  but  differs  by  0.4%  from 
Biele’s  value  which,  in  turn,  agrees  with  Tengelerfs  data  for  a  new 
40  mm  cannon. 

Based  on  the  values  cited  for  p  in  the  literature  and  on  this 
analysis,  we  can  safely  assume  that  the  uncertainty  in  the  term  Ap/p  is 
about  0.004.  It  will  produce  a  wavelength  error  on  the  order  of 
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AX  ,/A  ,  =  ±  0.00005  for  f.  =  15  GHz,  and 
gl  gl  1 


AX  0/A  _  =  ±  0.00233  for  f,,  =  2.725  GHz.  (6.2.10) 

g2  g2  2 


The  resulting  errors  in  projectile  travel  are  about  0.1  mm  and  4.7  mm 
for  one  meter  of  projectile  travel,  respectively,  accumulating  to  0.44  mm 
and  20.5  mm  over  the  total  in-bore  travel  for  the  105  mm  M68  tank  gun. 

This  estimate  shows  the  importance  of  knowing  the  exact  value  of  p,  if 
we  use  the  single  mode  approach.  With  today’s  instrumentation,  recording, 
and  data  analysis  techniques  there  is  no  excuse  for  not  determining  the 
electromagnetic  diameter  for  each  gun  system  in  the  inventory. 

6 . 3  Effect  of  Multiple  Propagating  Modes 

The  third  term  in  equation  (6.0.1),  JS^/(1  -  S^)](AR/R),  is  identi¬ 
cally  zero  for  the  single  mode  method.  For  the  multimode  technique, 
however,  we  must  consider  the  modulating  effect  of  the  higher  modes  on 
the  lowest  propagating  wave  mode  in  the  waveguide.  Since  most  of  the 
electromagnetic  field  energy  is  usually  partitioned  into  the  lowest 
propagating  wave  mode  as  indicated  by  the  amplitude  variation  of  the  shell 
push  data,  the  interference  with  the  higher  wave  modes  will  cause  the 
apparent  wavelength  to  fluctuate  between  two  bounds.  The  situation  may 
be  aggravated  because  of  the  geometry  and  angular  motion  of  the  surface 
reflecting  the  microwave.  Fbwever,  the  experiment  can  always  be  set  up 
such  that  the  energy  partition  requirement  is  fulfilled.  The  error  in 
the  displacement  measurement  will  be  limited  by  a  constant  throughout  the 
length  of  the  tube.  The  10  GHz  shell  push  data,  for  example,  exhibit  a 
maximum  accumulated  error  in  distance  which  is  for  a  flat  nosed  projectile 
less  than  13%  of  one  wavelength  at  any  instant  of  time.  For  our 
sensitivity  analysis,  we  can  safely  assume  that 


|  AX  /A  |  =  [S2/(l  -  S2)]  |  AR/R  |  <0.15  A  ,  (6.3.1) 

o  o  o 

limiting  the  uncertainty  in  projectile  travel  due  to  the  multimode  approach 
to  3  mm. 

By  using  a  data  analysis  technique  which  would  permit  decomposition 
of  the  recorded  time  series  into  individual  mode  components,  this 
oscillating  type  of  error  could  practically  be  kept  zero.  Hbwever,  no 
such  data  analysis  approach  is  yet  available. 

VII.  SUMMARY 

Combining  all  contributions  from  possible  errors  as  uncertainties 
in  the  parameters  of  the  guide  wavelength  equation,  we  can  establish 
a  bound  for  the  accumulative  error  in  the  measurement  of  projectile 
travel  via  microwave  interferometry  for  the  105  mm  M68  tank  gun  example, 
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TABLE  7.1.  UNCERTAINTIES  IN  THE  MEASUREMENT  OF  PROJECTILE  TRAVEL 


PARAMETERS  IN  THE  GUIDE  WAVELENGTH  EQUATION 

ESTIMATED  UNCERTAINTY  IN 
PROJECTILE  TRAVEL,  Imml 

15  GHz 

2.725  GHz 

FREQUENCY 

3.7 

3.7 

REFRACTIVE  INDEX: 

(a)  vacuum  approximation 

4.0 

6.7 

(b)  compressed  air  column 

2.0 

3.2 

CUTOFF  WAVELENGTH: 

(a)  bore  wear 

.  1 

4.9 

(b)  correction  for  rifling 

(c)  multimode  technique 

.5 

3.0 

2Q.5 

ACCUMULATIVE  ERROR  “ 

13.3 

39.0 

ACC.  ERROR/TOTAL  PROJ.  IN-BORE  TRAVEL 

.30% 

.89% 

These  errors  are  based  on  uncertainities  in  our  current  measurement  and 
data  analysis  procedures.  With  relative  moderate  efforts  we  can  reduce 
the  uncertainty  in  projectile  travel  at  least  by  one  order  of  magnitude: 

•  The  carrier  frequency  of  the  microwave  can  be  monitored  with  a 
frequency  counter  throughout  the  experiment.  This  would  practic¬ 
ally  eliminate  any  error  in  the  frequency. 

•  There  is  no  justification  in  using  the  vacuum  approximation  setting 
the  refractive  index  equal  to  one.  We  can  determine  the  index 

of  refraction  for  air  either  directly  by  measuring  it  with  a 
refractometer  or  indirectly  by  calculating  it  from  recorded 
atmospheric  parameters  such  as  pressure,  temperature  and  humidity. 

•  The  error  in  the  refractive  index  due  to  the  compressed  air  column 
with  and  without  blowby  gas  can,  in  principle,  be  determined 
experimentally.  However,  a  major  program  would  be  necessary  to 
accomplish  this,  especially  if  we  desire  quantitative  information 
on  obturation. 

•  The  effect  of  bore  wear  is  relatively  small  for  the  15  GHz  case 
and  can  be  neglected  for  all  practical  purposes.  For  the  single 
mode  measurement  approach,  however,  the  error  introduced  by 
bore  wear  may  not  be  negligible.  The  variation  of  bore  diameter 
as  function  of  tube  length  can  be  measured  by  star  gauging,  for 

example,  with  sufficient  accuracy  and  accounted  for  in  the’ data 
analysis . 
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•  The  uncertainty  in  the  electromagnetic  diameter  produces  an  error 
in  the  measurement  of  projectile  travel  which  is  relatively  small 
for  the  15  GHz  multimode  case  but  represents  the  largest  contri¬ 
bution  for  the  2.725  GHz  single  mode  case.  The  values  cited  in 
the  literature  for  it  are  too  far  apart.  Hence,  the  correction 
factor  for  the  helical  structure  of  the  bore  should  experimentally 
be  determined  for  each  gun  system  in  the  inventory  as  a  function 
of  the  microwave  carrier  frequency  with  an  accuracy  sufficient  to 
reduce  its  uncertainty  by  two  orders  of  magnitude.  I  t  is  a 
relatively  easy  task  to  design  instrumentation  which  uses  an 
infrared  laser  interferometer  as  a  displacement  reference  and 
readily  permits  the  required  measurement.  I  f  this  measurement 

is  done  as  an  integral  part  of  the  firing  experiment,  the 
correction  factor  given  as  function  of  axial  displacement  could 
even  include  the  effect  of  bore  wear. 

•  The  error  introduced  by  the  multimode  measurement  technique  can 
be  eliminated  by  going  to  the  single  mode  technique  or  can  be 
significantly  reduced  by  the  use  of  proper  data  analysis  techniques. 
Though  the  multimode  measurement  technique  requires  a  more 
complex  data  analysis  procedure  than  the  single  mode  method, 

it  has  the  intrinsic  advantage  that  it  allows  us  to  reduce  the 
errors  in  wavelength  introduced  by  parameter  uncertainties  or 
variations  by  sufficiently  increasing  the  microwave  carrier 
frequency. 

The  sensitivity  analysis  shows  that  the  currently  used  measurement 
procedures  lead  to  uncertainties  in  projectile  in-bore  travel  which 
exceed  the  less  than  one  millimeter  accuracy  requirement  for  the 
analysis  of  projectile-cannon  interface  problems.  However,  the  analysis 
also  exhibits  that  the  uncertainty  in  projectile  in-bore  travel  can  be 
reduced  to  less  than  one  millimeter,  if  we  employ  the  microwave 
measurement  technique  diligently  and  ensure  proper  determination  of 
all  parameters  in  the  guide  wavelength  equation. 


* 
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